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^ ■ Abstract 
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, Gutman et al. introduced the concepts of energy (^{G) and Laplacian energy 

Sl{G) for a simple graph G, and furthermore, they proposed a conjecture that for 
every graph G, ^{G) is not more than Si{G). Unfortunately, the conjecture turns 
out to be incorrect since Liu et al. and Stevanovic et al. constructed counterexam- 
Q . pies. However, So et al. verified the conjecture for bipartite graphs. In the present 

^ \ paper, we obtain, for a random graph, the lower and upper bounds of the Laplacian 

(-H \ energy, and show that the conjecture is true for almost all graphs. 
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\o 1 Introduction 

^ [ Throughout this paper, G denotes a simple graph of order n. The eigenvalues Ai, . . . , A„ of 
Q> I the adjacency matrix A{G) = {aij)nxn are said to be the eigenvalues of G. In chemistry, 
^ ' there is a closed relation between the molecular orbital energy levels of vr-electrons in 
conjugated hydrocarbons and the eigenvalues of the corresponding molecular graph. For 
^ ' the Hiichkel molecular orbital approximation, the total vr-electron energy in conjugated 

■ hydrocarbons is given by the sum of absolute values of the eigenvalues corresponding to 
the molecular graph G in which the maximum degree is not more than 4 in general. In 
1970s, Gutman [9] extended the concept of energy to all simple graphs G, and defined 
that 

n 

^(G) = ^|A.|, 

where Ai, . . . , A„ are the eigenvalues of G. Evidently, one can immediately get the energy 
of a graph by computing the eigenvalues of the graph. It is rather hard, however, to 
compute the eigenvalues for a large matrix, even for a large symmetric (0,l)-matrix like 
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A{G). So many researchers established a lot of lower and upper bounds to estimate 
the invariant for some classes of graphs. For further details, we refer readers to the 
comprehensive survey [H]. But there is a common flaw for those inequalities that only 
a few graphs attain the equalities of those bounds. Consequently we can hardly see 
the major behavior of the invariant (S'{G) for most graphs with respect to other graph 
parameters (|y((j')|, for instance). In the next section, however, we shall present an exact 
estimate of the energy for almost all graphs by Wigner's semi-circle law. 

In spectral graph theory, the matrix L(G) = D(G) — A{G) is called Laplacian matrix 
of G, where D(G) is a diagonal matrix in which da equals the degree dcivi) of the vertex 
Vi, i = 1, . . . ,n. Gutman et al. [I2] introduced a new matrix L(G') for a simple graph G, 
i.e., 

n n 

L(G) = L(G) - dG{vi)/nln = L(G) " 2 5Z 

i=l i=l j>i 

where I„ is the unit matrix of order n, and defined the Laplacian energy S'l{G) of G, i.e., 

n 

Sl{G) = Y,U. 
1=1 

where d, . . . ,(n are the eigenvalues of L(G). Obviously, we can easily evaluate the Lapla- 
cian energy S'l{G) if we could obtain the eigenvalues of L(G). In Section 3 we shall 
establish the lower and upper bounds of the Laplacian energy for almost all graphs by 
exploring the spectral distribution of the matrix L(G„(p)) for a random graph Gn{p) 
constructed from the classical Erdos-Renyi model (see [3]). 

In a recent paper [10], Gutman et al. proposed the following conjecture concerning 
the relation between the energy and the Laplacian energy of a graph. 

Conjecture 1. Let G be a simple graph. Then S'{G) < S'i{G). 

Unfortunately, the conjecture turns out to be incorrect. In fact, Liu et al. [H] and 
Stevanovic et al. [19] constructed two classes of graphs violating the assertion. However, 
So et al. [T7] proved that the conjecture is true for bipartite graphs. We shall show 
that the conjecture above is true for almost all graphs by comparing the energy with the 
Laplacian energy of a random graph in the third section. 

2 The energy of Gn(p) 

In this section, we shall formulate an exact estimate of the energy for almost all graphs 
by Wigner's semi-circle law. 

We start by recalling the Erdos-Renyi model Qn{p) (see [3]), which consists of all 
graphs with vertex set [n] = {1,2, . . . ,n} in which the edges are chosen independently 
with probability p = p{n). Apparently, the adjacency matrix A(G„(p)) of the random 
graph Gn{p) G Gn{p) is a random matrix, and thus one can readily evaluate the energy of 
Gn{p) once the spectral distribution of the random matrix A{Gn{p)) is known. 

In fact, the research on the spectral distributions of random matrices is rather abun- 
dant and active, which can be traced back to [21] . We refer readers to [H [6], [15] for an 
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overview and some spectacular progress in this field. One important achievement in that 
field is Wigner's semi-circle law which characterizes the limiting spectral distribution of 
the empirical spectral distribution of eigenvalues for a sort of random matrix. 

In order to characterize the statistical properties of the wave functions of quantum 
mechanical systems, Wigner in 1950s investigated the spectral distribution for a sort of 
random matrix, so-called Wigner matrix, 

X„ := (Xy), l<i,J<n, 

which satisfies the following properties: 

independent random variables with 

• the Xjj's have the same distribution Fi, while the Xi/s [i 7^ j) are to possess the 
same distribution F2; 

• "Var(xij) = < 00 for all 1 < i < j < n. 

We denote the eigenvalues of X„ by Ai,„, A2,n, • • • , ^n,n, and their empirical spectral dis- 
tribution (ESD) by 

(x) = #{Ai,„ I Xi^n <x, i = l,2,...,n}. 

Wigner [22l [23] considered the limiting spectral distribution (LSD) of X„, and obtained 
the semi-circle law. 

Theorem 1. Let X„ be a Wigner matrix. Then 

lim $„-i/2x„(a;) = $(x) a.s. 

n—>oo 

i.e., with probability 1, the ESD $n-i/2x„(^) converges weakly to a distribution $(x) as n 
tends to infinity, where $(a;) has the density 

Remark. It is interesting that the existence of the second moment of the off-diagonal 
entries is the necessary and sufficient condition for the semi-circle law, and there is no 
moment requirement on the diagonal elements. Furthermore, we can get more information 
about spectra of Wigner matrices. Set fit = J xdFi {i = 1, 2) and 

Xn = X„ — filln — /i2(Jr!. ~ In); 

where is the all I's matrix. One can easily check that each entry of X.„ has mean 0. 
By means of Wigner's trace method, one can show that the spectral radius p(n"^/2X„) 
converges to 2(72 with probability 1 as n tends to infinity (see Theorem 2 in [8], for 
instance). For further comments on Wigner's semi-circle law, we refer readers to the 
extraordinary survey by Bai [T]. 
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Following the book [3], we will say that almost every (a.e.) graph in Gn{p) has a certain 
property Q if the probability that a random graph Gn{p) has the property Q converges 
to 1 as n tends to infinity. Occasionally, we shall write almost all instead of almost every. 
It is easy to see that if Fi is a pointmass at 0, i.e., Fi{x) = 1 for x > and Fi{x) = 
for X < 0, and F2 is the Bernoulli distribution with mean p, then the Wigner matrix X„ 
coincides with the adjacency matrix of Gn{p)- Obviously, cr2 = \/p{^ ~ p) in this case. 

To establish the exact estimate of the energy £'{Gn{p)) for a.e. graph Gn{p), we 
first present some notions and assertions. In what follows, we shall use A to denote the 
adjacency matrix A(G„(p)) for convenience. Set 

A = A-p(J„-I„). 

Evidently, A is a Wigner matrix. By means of Theorem [1], we have 

lim $„-i/2A(a;) = ^{x) a.s. (1) 

n— >oo 

It is easy to check that each entry of A has mean 0. According to the remark above, 

lim p{n^^l^'X) = 2a2 a.s. (2) 

We further define the energy (f (M) of a matrix M as the sum of absolute values of 
the eigenvalues of M. By virtue of Equation ([T]) and ([2]), we shall formulate an estimate 
of the energy (S'{A), and then establish the exact estimate of (S'{A) = S'{Gn{p)) using 
Lemma [2l 

According to Equation ([2]), for any given e > 0, there exists an integer N such that 
with probability 1, for all n > the spectral radius p(n~^/^A) is not more than 2a2 + e. 
Since the density of is bounded on M, invoking Equation ([T]) and bounded 
convergence theorem yields that for all n > N, 

/i>2a2+e 
\x\d^j^-i/2x{x) = lim / \x\d^^-i/2-x{x) a.s. 

\x\d(^{x) a.s. (3) 

■2o-2-e 

\x\d^{x). 

We now turn to the estimate of the energy (S'{A). Suppose Ai, . . . , A„ and A^, . . . , A„ 
are the eigenvalues of A and n~^/^A, respectively. Clearly, X]r=i l-^^l ~ ''^^^^ SiLi 
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Equation we can deduce that 



1=1 



i=l 

j kM*n-l/2A(a^) 

j a.s. (n oo) 



2(72 



xU/4(To — (ix 



Therefore, with probabihty 1, the energy S" (A) enjoys the equation as follows: 

^ (A) = n^l' A ^^(13^ + 0(1)^ . 

We proceed to investigate (S'{A) = S'iGnij))) and present the following result due to 
Fan. 

Lemma 2 (Fan [7]). LetX, Y, Z he real symmetric matrices of order n such that^+Y = 
Z, then 

n n n 

J]|A,(X)| + ^|A,(Y)|>^|A,(Z)| 

i=l 1=1 i=l 

where Xi(M.) {i = 1, ■ ■ ■ ,n) is an eigenvalue of the matrix M. 

It is not difficult to verify that the eigenvalues of the matrix J„ — I„ are n — 1 and 
— 1 of ?i — 1 times. Consequently (S'{3n — In) = 2(n — 1). One can readily see that 
f?(p(J„-I„)) =pf?(J„-I„). Thus, 

^{p{3n-In)) =2p{n-l). 

Since A = A + p(J„ — it follows from Lemma [2] that with probability 1, 

^{A) < ^(A)+^(p(J„-I„)) 

Vp(1-p) + o(1) ) +2|>(r2-l). 



Consequently, 



lim (A)/n3/2 < —y^p{l-p) a.s. (4) 
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On the other hand, since A = A +p( — (J„ — !«)), we can deduce by Lemma [2] that with 
probabihty 1, 

<f(A) > ^(A) -^(p(-(J„-I„))) 
= ^(A) -^(p(J„-I„)) 

= ^3/2 -p) + o(1)^ - 2p{n - 1). 

Consequently, 



hm ^(A)/n3/2 > — -p) a.s. (5) 

n^oo 37r 

Combining equations (jl]) with ([5]), we have 

^(A)=n3/2 (J-^p(i-p) + o{l)^ a.s. 

RecaUing that A is the adjacency matrix of Gn{p), we thus obtain that a.e. random graph 
Gnip) enjoys the equation as follows: 



Note that for p = |, Nikiforov in [TB] got the above equation. Here, our result is for any 
probability p, which could be seen as a generalization of his result. 



3 The Laplacian energy of Gn(p) 

In this section, we shall establish the lower and upper bounds of the Laplacian energy of 
Gn{p) by employing the LSD of Markov matrix. Finally, we shall show that Conjecture 
[T] is true for almost all graphs by comparing the energy with the Laplacian energy of a 
random graph. 



3.1 The limiting spectral distribution 

We begin with another random matrix we are interested in. Define a random matrix 
M„ = X„ — D„ to be a Markov matrix if X„ is a Wigner matrix such that Fi is the 
pointmass at zero, and D„ is a diagonal matrix in which da = ^ = 1, . . . ,n. 

The matrix is introduced as the derivative of a transition matrix in a Markov process. 
Bryc et al. in [3] obtained the LSD of Markov matrix. Define the standard semi-circle 
distribution $0,1(2;) of zero mean and unit variance to be the measure on the real set of 
compact support with density 0o,i(3;) = 2^ V4 — x'^ l|a;|<2- 

Theorem 3 (Bryc et al. [S]). Let M„ he a markov matrix such that J xdF2{x) = and 
o"2 = 1. Then 

lim $„-i/2M„(a;) = ^(x) a.s. 
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where is the free convolution of the standard semi-circle distribution ^q i{x) and 

the standard normal measures. Moreover, this measure is a non-random symmetric 
probability measure with smooth bounded density, and does not depend on the distribution 
of the random variable Xij. 

Remark. To prove the theorem above, Bryc et al. employ the moment approach. In 
fact, they show that for each positive integer k, 



hm 

n— >oo 



j d^^-^/2^^{x) = j dm{x) a..s. (6) 



For two probabihty measures and there exists a unique probabihty measure /x ffl 
called the free convolution of /x and p. This concept introduced by Voiculescu [20] via 
C*-algebraic will be discussed in detail in the second part of this section. 

Let Gn{p) be a random graph of Gn{p)- Set a = \/p{l — p). One can easily see that cr^ 
is the variance of the random variable aij {i > j) in A(G„(p)). To state the main result 
of this part, we present a new matrix Li as follows: 

U = U{Gn{p)) = L{Gn{p))+p{3n-ln) (7) 

D(G„(p))-2^^a,,/nlJ - (A(a„(p)) - p(J„ - I„)) . 

i=l j>i / 

The following result is concerned with the LSD of Li. 
Theorem 4. Let Gn{p) be a random graph of Qn{p)- Then 

lim *(<xV7r)-iLi(a;) = ^(x) a.s. 

n— >cxD 

To prove the theorem above, we introduce an auxiliary matrix as follows: 

U = U{Gn{p)) = L(G„(p))-(r2-l)pI„ + p(J„,-I„) 

= (D(G„(p)) - (n - - (A(G„(j9))-p(J„-I„)). 

First of all, one can readily see that L2 is a Markov matrix in which the Wigner 
matrix is — A(G'„(p)) +p(Jn — In) and the diagonal matrix is — D(G'„(p)) + (n — l)pln- 
Furthermore, the off-diagonal entries of o"~^L2 have mean and variance 1. Since the 
LSD does not depend on the random variables Xjj, Theorem [3] yields 

lim $(^^)-iL2(a;) = ^(a;) a.s. 

n— >oo 

In what follows, we shall show that {a^/n)~^'Li and {a^/n)^^'L2 have the same 
LSD by which Theorem H] follows. To this end, we first estimate the difference 

[ay/n)~^ (Jji — L2) by Chernoff's inequality (see [13], pp. 26 for instance). 

Lemma 5 (Chernoff's Inequality). Let X be a random variable with binomial distri- 
bution Bi{n,p). Then, for any e > 0, 



P{\X-E{X)\ >e) <exp 



e2 



2{np - e/3) 
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Apparently, 



V 1=1 j>i ) 

Denote (o"\/^)~^(2 ^"^^^ X]j>i — {p- — by A„ for convenience. By means of 
Lemma 0, for any given e > 0, we have 



2^^aij/n- (n- l)p 



> e 



< 



P 



exp 



EE 

4 = 1 j>l 



> 



e ■ an 



3/2 



< exp 
= exp 



2(n(n - l)p + ean3/2/6) 



^(p + ecr/6) ■ ^2 



8(p + e(T/6) 

Therefore, by the first Borel-CanteUi lemma (see [2], pp. 59 for instance), we can deduce 



0" 



i=\ j>i 



a.s. (n oo) 



Furthermore, it is easy to see that A is an eigenvalue of {o^^/n) ^Li if and only if A + A„ 
is an eigenvalue of {ay/n)^^L2. By the definition of the ESD, it follows that 

= "^{aV^r^ui^ + An). (8) 

Clearly, for any e > 0, there exists such that |A„| < e a.s. for all n > N. Noting that 
'^(crv^)-iL2(^) is an increasing function, for all n> N, we have 

^{aV^r^hA^ - e) < $(<xV^)-iL2(a^ + A„) < ^(aV^)-iL2ix + c) a.s. 
Consequently, 

^{x-e) = lim $(^^). 112(3^ -e) 

< lim $(a^)-iL2(a; + A„) 

n^oo 

< lim <l>(^ /^)-iL2(x + e) = ^(x + e) a.s. 

n— >oo 

Moreover, since the density of \Ef(x) is smooth bounded, \E'(x) is continuous. Together 
with the fact that e is arbitrary, we conclude 



lim $(^^)-iLi(a;) = lim $ 



'x + A„) = ^(x) a.s. 



which completes the proof of Theorem |H 
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3.2 The bounds of (^L{Gn{p)) 

In this part, we shall establish the lower and upper bounds of S'L{Gn{p)) by employing 
Theorem m and then show that Conjecture [T] is true for almost all graphs at last. 

Let X be a random variable with the distribution '^{x). We start with an estimate of 
E |X| = j \x\d'^{x). Since '^{x) is the free convolution of the standard semi-circle distri- 
bution $0,1 (^) and the standard normal measure, let us investigate the free convolution 
in depth. Here, we follow the notation given by Voiculescu [21]. The Cauchy-Stieltjes 
transform of a probability measure is 



which is analytic on the complex upper half plane. For some a,(3 > 0, there exists a 
domain -Da,/? = {u + iv \ \u\ < av,v > (3} on which is univalent. For the image 
Gfj_{Da,i3), we can define the inverse function of in the area Ta^b = {u + iv \ \u\ < 
—av, —b<v<0}. And let R^{z) = K^{z) — Ij z. Then for probability measures /i and 
there exists a unique probability measure, denoted by ffl i/, on such that 

RiMv = Rfi + Ri/- 

The measure fiS u is said to be the free convolution of fi and u. 

In the above definition, the Cauchy-Stieltjes transform and inverse function may be 
difficult to compute in practice. Consequently, we do not compute E |X| directly. In 
what follows, we employ another definition of free convolution via combinatorial way (see 
[U [18]) applicable only to probability measures with all moments. 

For probability measure ^, set rrik = J x^fi{dx) and 

oo 

M^{z) = l + ^mfc2^ 

k=l 

Define a formal power series 

oo 

T,{z)=J2ckz'-' 

k=l 

such that 

M^{z) = l + zM^{z)T^{zM^{z)). 
Then, the free convolution of /i, u is the probability measure fiS u satisfying 

T^muiz) =T^{z) + T,{z). (9) 



It is not difficult to see that this definition is coincident with the analytical one (see [H 

Next, we calculate E|X| by the following result due to Bryc Let M^,„ = M^{z) 
mod z"'~^^, T^^n{,z) = T^{z) mod z"-^^ be the n-th truncations, i.e., M^,„ = 1 + ^^=1 rrikZ^ 
and T^,niz) = Y^ltl CkZ^~^. 
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Lemma 6 (Bryc ||3]). With M^^^z) = 1 and ci = ^^(0), we have 

M^^niz) = 1 + zM^,n-l(2)T^,n-l(2:M^,n-l(^)) Hiod n > 1, 

and 

1 1 d' 

Ck = 



k-lk\ dz^ M^--[z] 



2 = 

Therefore, combining with the formula ([9]), we can calculate the moments of yU ffl 
by the moments of v in recurrence. It is not difficult to verify that EX^ = 2 and 
EX^ = 9 (see [1] for details). Employing Cauchy-Schwartz inequality 



E(Xr)p < EX^ -EF 



we have 
and 

Therefore, 



EIXI < VEX2 



[¥jX'^f < E |X| ■ E < E |X| ■ VEX2 -EXl 
^^<E|X|<v^. 



3 

In what follows, we shall establish the lower and upper bounds of ^L(G'n(p)) by 
employing an estimate of the energy S{Li). We first investigate the convergence of 
/ |x|(i$(o-^)~iLj(x). Let I be the interval [—1,1]. By Theorem H] and the bounded 
convergence theorem, one can easily see that 

lim / \x\d^i„ Ax) = I \x\d'^{x) cL.s. (10) 

n^ooj^ J I 

We proceed to prove that 

lim / |z|c?$(^^)-iL^(x) = / \x\d-^{x) 

where J'^ = M \ /. Since cr~^ L2 is the Markov matrix such that the off-diagonal entries 
have mean and variance 1, we can deduce, by Equation ([6]), that 

lim / /;^wiL (x) = / x^(i\E'(x) a.s. (11) 

n^ooj ^ ^ ^ J 

According to the relation ([8]), we have 

= y"(x- A„)2d$(,^)-iL2(a;) 

= j x^ci$(^^)-iL2(x) - 2A„ j xd$(^07)-iL2(a;) 

+AlJ d<l>(,^)-iL2(x). 
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Since lim„^oo^n = a.s., Equation ffTTj) implies that 

lim / a;2rf$(,^)-iL^(x) = lim / x^d^(^^^yi^^{x) = / x^d-^{x) a.s. (12) 



Consequently, 



a.s. 

Lemma 7 (Billingsley [2] pp. 219). Let fi be a measure. Suppose that functions 
Ctrl) i>n, fn converges almost everywhere to functions a, b, f , respectively, and that an < fn < 
bn almost everywhere. If J Qndfi —yjadjji and J bndfi J bdjji, then J fndfi — ^ / fdfi. 

Suppose that 0(<jv^)-i Li(3^) is the density of $(crv^)-iLi(3^)- By virtue of Theorem 
m and Lemma [71 we can deduce by setting a„(x) = 0, 6„(x) = a;^0(o.^)-iL^(x) and 
fn{x) = |a;|0(^^)-iLi(x) that 

lim / |x|(i$(o.^)-iL^(x) = / |x|(i\l/(x) a.s. 

Combining the above equation with Equation (fTUj) . we have 

lirn^ j |x|(i$(^^)-iL^(a;) = j \x\d^{x) a.s. 

We are now ready to present an estimate of the energy (f (Li). By an argument similar 
to evaluate the energy ^(A), we have 

S{U)/an^l^ = J |x|d<l>(,^)-.L,(x) 

|a;|(i\&(a;) a.s. (n — > oo). 



Since 2^2/3 < E|X| < ^2, 
Consequently, 



o _ Q/9 ^ V - a.s. (n oo] 



+ o(l) j n""^ < ^(Li) < [V2a + o(l)) n^'^ a.s. (13) 

Employing the equation above, we can establish the lower and upper bounds of 
<S'L{Gn{p))- Note that S'liGnip)) = <^(L) according to the definition of the energy of 
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a matrix. So we turn our attention to the bounds of (S'(L). By means of Equation ([7]), 
we have 

Li = L +p{3n - In) and L = Li +p(I„ - J„). 
Thus, Lemma [2] yields that 

<f (Li) - ^(p(J„ - I„)) < ^(L) < ^(Li) + - J„)). 

Recalhng the fact that S'{p{3n — !«)) = <^(p(In — Jn)) = 2p(n — 1), Equation ( fT3l) imphes 
that 

n — 1) a.s. 

Therefore, we obtain the the lower and upper bounds of the Laplacian energy for almost 
all graphs. 

Theorem 8. Almost every random graph Gn{p) satisfies 



+ 0(1) ■ < ^i(G„(p)) < (y2a + o(l)) 



Since a.e. random graph Gn{p) satisfies 

^{G^jp)) 8 2v^ ^L(gn(p)) 
lim ^75 — = —a < -^(T < hm — , 

we thus establish the result below. 

Theorem 9. For almost every random graph Gn{p), S{Gn{p)) < (S'i{Gn{p))- 
By virtue of the theorem above. Conjecture [T] is true for almost all graphs. 

Acknowledgement: The authors are very grateful to the referees for detailed suggestions 
and comments, which helped to improve the presentation of the manuscript significantly. 
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